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Abstract 



Current neutrino oscillation data indicate that 9 13 is not strongly suppressed 
and #23 might have an appreciable deviation from 7r/4, implying that the 
3x3 neutrino mixing matrix V does not have an exact \x-t permutation 
symmetry. We make a further study of the effect of \i-t symmetry breaking on 
the democratic flavor distribution of ultrahigh-energy (UHE) cosmic neutrinos 
at a neutrino telescope, and find that it is characterized by IV^J 2 — \V ri \ 2 which 
would vanish if either 6 23 = vr/4 and 6 13 = or 6 23 = 7r/4 and 5 = ±tt/2 held. 
We observe that the second-order fi-r symmetry breaking term A may be 
numerically comparable with or even larger than the first-order term A in the 
flux ratios : <j>J : <ff ~ (1-2A) : (1 + A + A) : (1 + A-A), if sin((9 23 -vr/4) 
and cos 5 have the same sign. The detection of the UHE V e flux via the 
Glashow-resonance channel V e e — > W~ — >• anything is also discussed by taking 
account of the first- and second-order \x-r symmetry breaking effects. 
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I. INTRODUCTION 



Current experimental data have convinced us that three known neutrinos can oscillate 
from one flavor to another [1], implying the existence of a mismatch between their flavor 
and mass eigenstates. Hence the lepton flavors must mix as the quark flavors, and this 
phenomenon can be described by using an effective 3x3 unitary matrix V. In the basis 
where the flavor and mass eigenstates of three charged leptons are identical, V provides a 
unique link between the neutrino flavor (u e , u T ) and mass {v x , u 2 , u 3 ) eigenstates: 
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(v el v e2 v e3 
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One may parametrize V in terms of three mixing angles and three phase angles as follows: 
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P v , (2) 



where c^- = cos9 ij: = sin^- (for ij = 12, 13,23), and P v = Diag{e* p , e ta , 1} is physically 
relevant only if massive neutrinos are the Majorana particles. It has been pointed out that 
|V^j| = \V Ti \ (for % = 1,2,3) holds exactly, if either the conditions 9 13 = and 6 23 = 7r/4 
or the conditions 5 = ±n/2 and 9 23 = 7r/4 are satisfied [2]. However, this interesting /i-t 
permutation symmetry must be broken: on the one hand, the fact that 9 13 is not strongly 
suppressed has recently been established in the Daya Bay and RENO reactor antineutrino 
oscillation experiments [3] ; on the other hand, the latest global fit of all the available neutrino 
oscillation data hints that the value of 9 23 might have an appreciable deviation from 7r/4 [4]. 
An important task in today's experimental neutrino physics is therefore to determine the 
strength of /i-t symmetry breaking, so as to help the theorists try different flavor symmetries 
and deeply understand the leptonic flavor mixing structure [5]. 

The neutrino telescopes (e.g., the running IceCube detector at the South Pole [6] and 
the proposed KM3NeT detector in the Mediterranean Sea [7]), which aim to observe the 
ultrahigh-energy (UHE) cosmic neutrinos and their flavor distributions, can serve as a novel 
probe of the /i-r symmetry breaking effects. It is well known that neutrino oscillations may 
map 0g : 0^ : 0^ = 1 : 2 : 0, the initial flavor ratios of UHE cosmic neutrino fluxes produced 
from p7 or pp collisions at a distant astrophysical source, into 0j : 0^ : 0^ = 1 : 1 : 1 at 
a neutrino telescope on the Earth [8] if there is the exact /i-r permutation symmetry. It is 
also known that such a democratic flavor distribution can be broken down to [9] 

0j:0j:0^(l-2A):(l + A):(l + A) , (3) 

where 

A = - sin 2 29 12 sin e — - sin 49 12 sin 9 13 cos 5 , (4) 

and e = 9 23 — 7r/4 for the parametrization of V given in Eq. (2). Hence A signifies the 
combined effects of /i-r symmetry breaking, and whether its magnitude can reach the 10% 
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level or not depends crucially on the sizes and signs of sin e and cos 5. Note that the analytical 
approximation made in Eq. (3) does not reflect the difference between 0^ and <f>^, which 
should be given by the terms proportional to sin 2 e, sin 2 9 13 and sine sin 9 l3 [10]. 

The purpose of this paper is to make a further study of /i-t symmetry breaking relevant to 
the detection of UHE cosmic neutrinos originating from a certain cosmic accelerator. First, 
we calculate the flavor distribution of UHE cosmic neutrinos at a neutrino telescope in a 
parametrization-independent way. Such an exercise allows us to generalize the approximate 
result in Eq. (3) to the following exact one: 



jT 0o 
a = y 



i + £K/(l^l 2 -IKJ 2 ) 



(5) 



for a — e, /j, and r, where O = 0^ + 0^ + 0^ is the total flux of UHE cosmic neutrinos 
and antineutrinos of all three flavors. Now it becomes transparent that a deviation from 
the democratic flavor distribution 0j : 0^ : 0j = 1 : 1 : 1 at a neutrino telescope is 
characterized by the fi-r symmetry breaking terms |V^| 2 — |U ri | 2 . In particular, we find 
that the difference between 0^ and 0^ is measured by (|V^| 2 — |U T j| 2 ) 2 . We also find that 
the second-order /i-t symmetry breaking term A, which is a function of sin 2 e, sin 2 9 l3 and 
sine sin 9 13 , may be numerically comparable with or even larger than A in the realistic flux 
ratios 0j : 0j : 0j ~ (1 - 2A) : (1 + A + A) : (1 + A - A) if sine and cos 5 have the 
same sign. Second, we reexamine the effect of fi-r symmetry breaking on the V e flux of 
E v ~ 6.3 PeV which is detectable via the well-known Glashow-resonance (GR) channel 
V e e — > W~ — > anything [11] at neutrino telescopes. Different from the previous result 
obtained Ref. [9], our new result for the GR-mediated V e events will include the second- 
order ji-r symmetry breaking effect. 



II. EFFECTS OF /i-t SYMMETRY BREAKING 

Let us define 0^ = 0^ + 0§ and 0^ = 0^ + 0^ (for a — e, /j,, r) throughout this paper, 
where 0^ (or 0j ) and 0| (or 0£ ) denote the v a and V a fluxes at a distant astrophysical 
source (or at a neutrino telescope), respectively. For most of the currently-envisaged sources 
of UHE cosmic neutrinos [12], a general expectation is that the initial neutrino fluxes are 
produced via the decay chain of charged pions and muons created from pp or p7 collisions 
and their flavor content can be expressed as 



{0e S 



S = {\, |, O^0 O , (6) 



where (f)f = 05 = 0| = 0, and O = 0^ + 0^ + 0^ is the total flux of neutrinos and 
antineutrinos of all flavors. Thanks to neutrino oscillations, the flavor distribution of such 
UHE cosmic neutrinos at a neutrino telescope is described by 

€ = <Pl + 4>l = E k,Piy P -+ v a) + -+ v a) \ = EE (\vj 2 \v,f4) . ( 7 ) 

where we have used 
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(8) 



Since the Galactic distance that the UHE cosmic neutrinos travel far exceeds the observed 
neutrino oscillation lengths, P{vq — > f a ) and P(pp — > v a ) are actually averaged over many 
oscillations and thus become energy- and distance-independent. Combining Eq. (7) with 
Eq. (6) and using the unitarity conditions of V, we find 



00 , x = 00 ,y |2 + |2_ | y |2\ 

o / / I at I \l eil 1 I v fill J o / , I v m I ^ 1 I v I l v ril y 



(9) 



This is a simple proof of Eq. (5). Of course, the relationship 0j + 0^ + 0^ = O holds. To 
be more explicit, we have 



i + EIK*I 2 (I^I 2 -IK* 

i 

i+ei^i 2 (i^i 2 -ikj 2 ) 

i 

i + ElKJ 2 (l^l 2 -lK*l 2 ) 



(10) 



from which we obtain the difference between 0^ and 0^ as 

«-^ = ^E(l^l 2 -IKJ 2 ) 2 



(11) 



Now it becomes quite transparent that the deviation of 0^ (for a = e,fj,,r) from o /3 is 
measured by |V^| 2 — IV^J 2 , and the difference between 0^ and 0^ is purely an effect governed 
by (|V^J 2 — l^-J 2 ) 2 - This exact and parametrization-independent observation is therefore 
useful for us to probe the leptonic flavor mixing structure via the detection of UHE cosmic 
neutrinos at neutrino telescopes. 

Considering that the \i-r symmetry of V is possible to be partly or softly broken, we define 
the following three fi-r symmetry breaking quantities and express them in terms of three 
neutrino mixing angles and the Dirac CP-violating phase in the standard parametrization 
of V as given in Eq. (2): 



- (sin 2 9 12 — cos 2 9 12 sin 2 13 ) cos 29 23 + sin 29 12 sin 29 23 sin 9 13 cos 5 , 
A 2 = | V^ 2 1 2 — | V t2 1 2 = (cos 2 9 12 — sin 2 9 12 sin 2 13 ) cos 29 23 — sin 29 12 sin 26> 23 sin 9 13 cos 5 , 



\ = 


W A 2 - 


IKi 


A 2 = 


W I 2 - 

1 >2l 


IK2 


A 3 = 


K,\ 2 - 


IK3 



= — cos 9 13 cos 26*. 



'23 



(12) 



We see that the relationship A 1 + A 2 + A 3 = holds exactly, as guaranteed by the unitarity 
of V. There are three special but interesting cases, in which the fi-r permutation symmetry 
is not completely broken: 
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• Case (A): 9 23 — > ir/4. In this limit, A i (for i — 1, 2, 3) can be simplified to 

A 1 = — A 2 = sin 29 12 sin # 13 cos 5 , 

A 3 = . (13) 

• Case (B): 9 13 — > 0. In this limit, the results in Eq. (12) are simplified to 

A 1 = tan 2 9 12 A 2 = sin 2 6* 12 cos 29 23 , 

A 3 = -cos 2^3 . (14) 

• Case (C): 5 — > ±ir/2. In this limit, we simply obtain 

A x = (sin 2 9 12 — cos 2 9 12 sin 2 9 l3 ^j cos 2# 23 , 
A 2 = (cos 2 9 12 — sin 2 9 12 sin 2 6> 13 ) cos2# 23 , 
A 3 = — cos 2 9 13 cos 29 23 . 



(15) 



Current experimental data tell us that case (A) is probably not true, case (B) is definitely not 
true, and case (C) remains an interesting possibility. In any case, we may make analytical 
approximations for 0^ (for a — e, /i, r) up to the accuracy of sin 2 e, sin 2 9 13 and sine sin 9 13 . 
Then a result similar to Eq. (3) is 

<fl : ft : f T ~ (1 - 2A) : (1 + A + A) : (1 + A - A) , (16) 

where A has been given in Eq. (4), and A is defined as 1 

A = (4 — sin 2 29 12 ^j sin 2 e + sin 2 26> 12 sin 2 9 13 cos 2 5 + sin 4# 12 sin e sin 9 13 cos 5 . (17) 

It is easy to show that A > holds for arbitrary values of 5, because the above expression 
can simply be transformed into A = 3sin 2 e+ (cos 29 l2 sine + sin 29 12 sin 9 13 cos 5) 2 , which 
may vanish if both sin e = and sin 9 13 = (or cos 5 = 0) hold. 

One may define three working observables at neutrino telescopes [13] and link them to 
the [i-t symmetry breaking quantities: 



R, 



R = 



0! 


i + E 




2 A, 
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V ■ 

ei 


2 A, 




i + E 


v.. 


2 A, 




2-E 
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V ■ 

ill 


2 A, 




i + E 
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V ri 


2 A, 




2-E 


V ri 


2 A, 



1 

2 

1 

2 



i + ^EIK,I 2 a, 



2 2 



'i + |i;ivy 2 A,W + §(A + A) , 



4(i + |i:ik,i 2 a,)4 + |(a-a). 



(18) 



1 This second-order perturbation term is consistent with the one obtained in Ref. [10], where the 
definition of e has an opposite sign. 
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Comparing Eq. (11) and Eqs. (16) — (18), we immediately obtain 




Hence A signifies the second-order effect of fi-r symmetry breaking, and the departure of 
R a (for a = e, //, r) from 1/2 is a clear measure of the overall \x-r symmetry breaking effects. 

To illustrate the possible size of \i-r symmetry breaking, we estimate A, (for % = 1, 2, 3), 
A and A by using the latest values of 9 12 , 9 13 , 9 23 and 5 obtained from a global analysis 
of the presently available neutrino oscillation data done by Fogli et al [4]. Our numerical 
results are listed in TABLE I, where both normal and inverted neutrino mass hierarchies 
are taken into account in accordance with Ref. [4]. Three comments are in order. 

• | AJ ~ 0.1 (for i = l,2, 3) holds when the best-fit values of 9 12 , 9 13 , 9 23 and 5 are used. 
Given the 2a intervals of the four input parameters, |AJ ~ is allowed except for A 3 
in the case of a normal neutrino mass hierarchy. It is therefore desirable to determine 
the departure of 9 23 from 7r/4 and that of 5 from ±n/2 in the future experiments. 

• Since the best-fit values of 9 23 and 5 lie in the ranges < 9 23 < ir/4 and n/2 < 5 < n 
respectively, sine and cos<5 have the same sign and thus the two terms of A in Eq. (3) 
significantly cancel each other. This significant cancellation leads to |A| < A at the 
percent level, implying a fairly good flavor democracy for 0j, 0^ and (fij. 

• When the 2a intervals of 9 12 , 9 13 , 9 23 and 5 are taken into account, the lower bound of 
A and the upper bound of A are about 0.1 in magnitude. Hence the flavor democracy 
of 0j, 0^ and 0^ can maximally be broken at the same level. Needless to say, an 
appreciable departure of R a (for a = e,/i,T) from 1/2 requires an appreciable effect 
of \x-r symmetry breaking. 

The rapid development of reactor antineutrino oscillation experiments implies that 9 l3 will 
soon become the best known angle of V. In this case, improving the precision of 9 23 and 
determining the value of 5 turn out to be two burning issues which will allow us to pin down 
the leptonic flavor mixing structure including CP violation. 

III. ON THE GLASHOW RESONANCE 

Let us proceed to look at the effect of /i-t symmetry breaking on the flavor distribution 
of UHE cosmic neutrinos at a neutrino telescope by detecting the V e flux from a distant 
astrophysical source through the GR channel V e e — > W~ — > anything [11]. This reaction can 
take place over a narrow energy interval around the V e energy E^ R My / /2m e w 6.3 PeV, 
and its cross section is about two orders of magnitude larger than the cross sections of V e N 
interactions of the same V e energy [14]. A measurement of the GR reaction is important 
in neutrino astronomy because it may serve as a sensitive discriminator of UHE cosmic 
neutrinos originating from and pp collisions [15-18]. We hope that a neutrino telescope 
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may measure both the GR-mediated V e events and the v + V events of charged-current 
interactions in the vicinity of E^ R , and their ratio can be related to the ratio of the V e flux 
to the v.. and V„ fluxes entering the detector: 



T 



Here we follow Ref. [9] to reexamine the fx-r symmetry breaking effect on -R GR . 

The initial UHE cosmic neutrino fluxes are produced via the decay chain of charged pions 
and muons created from pp or pj collisions at a cosmic accelerator, and thus their flavor 
distribution can be expressed as 

\ W ' \ ' \ ' \ ' ° ' °l ^° ( pp coiiisi ° ns ) > 

[ 1 3 , , - , - , , j O (pry collisions) . 

In either case the sum of 0^ and 0§ (for a = e,/j,,r) is consistent with c/r in Eq. (6). 
Thanks to neutrino oscillations, the V e flux at a neutrino telescope is given by 

<f>l = E [4p(^ -+ v e )] = EE (KmM,) • (22) 

To be explicit, 

0j e (H = ^(l + ElKJ 2 A^^(l-2A) , 

<£>7) = f E \VJ 2 K t \ 2 - I [sin 2 2£ 12 - 4A + 2 (l + cos 2 29 12 ) sin 2 13 ] . (23) 

i 

Since the expression of 0^ can be found in Eq. (10), it is straightforward to calculate R GR 
by using Eqs. (20) and (23) for two different astrophysical sources: 

r, / X 1 3 . 1 — 

Rgr(pp) - 2 ~ 2 2 ' 

^ GR (P7) * ^ - i^^A - ^ A + sm 2 * 13 . (24) 

We see that the deviation of R GR (pp) from 1/2 and that of Rq R (p~y) f rom sin 2 2# 12 /4 are 
both controlled by the effects of fi-r symmetry breaking, which can maximally be of 0(0.1). 
As discussed in Refs. [15] and [17], the IceCube detector running at the South Pole has a 
good discovery potential to measure Rq R (pp) after several years of data accumulation. In 
comparison, it seems more difficult to probe the GR-mediated UHE u e events originating 
from the pure p7 collisions at a cosmic accelerator. 

Of course, one may also consider some other possible astrophysical sources of UHE cosmic 
neutrinos, such as the neutron beam source [19] with {0^ : 0^ : 0^} = {1:0:0} and the 
muon-damped source [20] with {0^ : 0^ : 0^} = {0:1:0}, to study their flavor distributions 
at neutrino telescopes and probe the effects of /i-t symmetry breaking. 
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IV. SUMMARY AND FURTHER DISCUSSIONS 



With the development of several neutrino telescope experiments, a lot of interest has 
recently been paid to the flavor issues of UHE cosmic neutrino fluxes and whether they are 
detectable in the foreseeable future [21]. In view of the fact that the smallest neutrino mixing 
angle 8 13 is not strongly suppressed [3] and the hint that the largest neutrino mixing angle 
#23 might have an appreciable departure from 7r/4 [4], we have carried out a further study 
of the effect of fi-r symmetry breaking on the presumably democratic flavor distribution 
of UHE cosmic neutrinos at a neutrino telescope. Our new results are different from the 
previous ones in the following aspects: 

• An exact and parametrization-independent expression for the fluxes of UHE cosmic 
neutrinos at a neutrino telescope is obtained in Eq. (5) or Eq. (10), and it clearly 
shows the /i-r symmetry breaking effect measured by \V^\ 2 — \V ri \ 2 (for i = 1,2,3). 
In addition, the difference between 0^ and 0!jF is the pure second-order /i-r symmetry 
breaking effect proportional to (|V^| 2 — |U Ti | 2 ) 2 . 

• The first- and second-order /x-r symmetry breaking effects, characterized respectively 
by A and A in the flux ratios 0j : 0^ : <pj ~ (1 - 2A) : (1 + A + A) : (1 + A - A), may 
be numerically at the same order of magnitude. This observation is particularly true 
when sine and cos<5 have the same sign such that the two first-order jj,-r symmetry 
breaking terms of A cancel each other and then lead us to |A| < A. 

• A > holds for arbitrary values of 5, and its contributions to the flux ratios R GR (pp) 
and i? GR (p7) for the GR-mediated V e events are taken into account in Eq. (24). The 
term proportional to sin 2 9 13 is also included in the expression of -R GR (j97), but such a 
term does not appear in R GR (pp). 

Because of the poor data on e and 5, whether the magnitude of A and A (or one of them) 
can be as large as about 10% remains an open question. It is also possible that both of them 
are at the 1% level, as shown in TABLE I with the present best-fit results of 6 23 and 5 [4]. If 
|A| and A are finally confirmed to be really small, then an approximately democratic flavor 
distribution 0j : 0^ : 0^ ~ 1 : 1 : 1 will show up at neutrino telescopes for the UHE cosmic 
neutrinos originating from pp and (or) pj collisions at a distant astrophysical source. 

So far we have omitted the uncertainties associated with the initial neutrino fluxes at 
a given astrophysical source. A more careful analysis of the flavor ratios of UHE cosmic 
neutrinos originating from pp or pj collisions yields 0^ : 0^ : 0^ ~ 1 : 2 (1 — rj) : with 
i] ~ 0.08 [22]. If this uncertainty is taken into account, Eq. (5) or Eq. (9) will change to 




EI^PKl^l 2 -!^! 2 )- 2 ^^! 2 ] (25) 



for a — e, fx and r. This result clearly shows that the 77-induced correction is in general 
comparable with (or even larger than) the effect of /i-r symmetry breaking. On the other 
hand, there exist the uncertainties associated with the identification of different flavors at 
the neutrino telescope. Taking the IceCube detector as an example, Beacom et al [23] have 
pointed out that the relative experimental error for determining the ratio of the muon track 
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to the non-muon shower is typically £ ~ 20%, depending on the event numbers. This 
estimate means that the ratio R in Eq. (18) may in practice be contaminated by £, which 
is very likely to overwhelm the [i-r symmetry breaking effect (A + A). 

We also admit that our treatment does not take account of the other complexities and 
uncertainties associated with the origin of UHE cosmic neutrinos, such as their energy 
dependence, the effect of magnetic fields and even possible new physics [24]. It is still too 
early to say that we have correctly understood the production mechanism of UHE cosmic rays 
and neutrinos from a given cosmic accelerator. But the progress made in the measurement 
of lepton flavor mixing parameters is quite encouraging, and it may finally allow us to well 
control the error bars from particle physics (e.g., the effect of /i-r symmetry breaking) and 
thus concentrate on the unknowns from astrophysics (e.g., the initial flavor composition of 
UHE cosmic neutrinos). We believe that any constraint on the flavor distribution of UHE 
cosmic neutrinos to be achieved from a neutrino telescope will be greatly useful in diagnosing 
the astrophysical sources and in understanding the properties of neutrinos themselves. Much 
more efforts are therefore needed to make in this direction. 
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TABLES 



TABLE I. Possible sizes of fi-T symmetry breaking at neutrino telescopes. The values of # 12) 
# 13 , #23 and 5 are taken from the latest global analysis of currently available neutrino oscillation data 
done by Fogli et al [4], where both normal and inverted neutrino mass hierarchies are considered. 



Normal hierarchy 


Best fit 


2d rang 


;e 




sin 2 6> 12 


3.07 x KT 1 


(2.75 ■ ■ 


■ 3.42) x 10" 


-i 


sin 2 6> 13 


2.45 x KT 2 


(1.81 •• 


■3.11) x 10- 


-2 


sin 2 # 23 


3.98 x KT 1 


(3.50 ■■ 


■ 4.75) x 10- 


-1 


5 


0.89 X 7T 


• • • 2vr 






Ai 


-7.38 x 10" 2 


-0.151 


••• + 0.256 




A 2 


+2.73 x lO" 1 


-0.135 


••• + 0.365 




A 3 


-1.99 x lO" 1 


-0.295 


0.048 




A 


-1.74 x 10~ 2 


-0.096 


• • • + 0.026 




A 


+6.23 x 10~ 2 


- 0.120 




Inverted hierarchy 


Best fit 


2d rang 


;e 




sin 2 6> 12 


3.07 x 10- 1 


(2.75 ■ ■ 


■ 3.42) x 10" 


-1 


sin 2 6> 13 


2.46 x 10" 2 


(1.83 •• 


•3.13) x 10" 


-2 


sin 2 (9 23 


4.08 x lO" 1 


(3.55 •• 


• 6.27) x 10" 


-1 


5 


0.90 X 7T 


■ ■ ■ 2ir 






Ai 


-8.19 x 10~ 2 


-0.244 


• • • + 0.254 




A 2 


+2.61 x 10" 1 


-0.335 


• • • + 0.359 




A 3 


-1.79 x lO" 1 


-0.285 


• • • + 0.249 




A 


-1.28 x 10~ 2 


-0.094 


• • • + 0.087 




A 


+5.56 x 10- 2 


0--- 0.115 
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